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STABILITY  IN  NEUTRAL  EQUATIONS 
by 


Jack  K.  Hale  and  Pedro  Mar tine 2-Amores 


Abstract 


Coupled  systems  of  differential-difference  and  ordinary 
difference  equations  occur  in  various  applications  including  the 
theory  of  transmission  lines  [1]  and  gas  dynamics  [2] . Stability 
of  linear  systems  has  been  discussed  by  Brayton  [1]  using  Laplace 
transform  and  the  problem  of  absolute  stability  by  Rasvan  [12] 
using  the  frequency  domain  method  of  Popov. 

In  this  paper,  the  same  problems  are  discussed  by  the 
following  method.  By  differentiating  the  difference  equation, 
one  obtains  a system  of  neutral  differential-difference  equations. 
The  desired  solutions  of  the  original  problem  are  obtained  by 
restricting  the  initial  data  to  lie  on  certain  manifolds  in  the 
space  of  all  initial  data.  In  this  way,  this  class  of  problems 
can  be  treated  in  a natural  manner  using  the  known  methods  of 
neutral  equations.  Generalizations  to  arbitrary  functional 
differential  equations  is  also  immediate  when  this  approach  is 
employed . 
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(1.2) 


at  D(xt)  = L(xt} 


A solution  x = x($)  through  (j>  e C at  t=0  is  a 
continuous  function  taking  [-r,A) , A > 0,  into  Rn  such  that 
Xq  = <j>,  D(xt)  is  continuously  differentiable  on  [0,A)  and 
equation  (1.2)  is  satisfied  on  this  interval.  It  follows  from 
Hale  and  Mayer  [9]  that  a solution  through  <j>  exists  on  I-r,°°) , 
is  unique  and  depends  continuously  in  (J>. 

If  T(t):  C C,  t > 0,  is  defined  by  T(t)<{>  = xt(<}))/ 
then  T(t),  t > 0 is  a strongly  continuous  semigroup  with  in- 
finitesimal generator  A:  ^(A)  C,  A<J>(0)  = $(6),  -r  < 0 < 0, 


Q>{&)  = e C:  <f>  e C,  D$  = L<p). 


The  spectrum  cr(A)r  of  A consists  of  all  X which 


satisfy  the  characteristic  equation 


Cl. 3)  det  A (X)  = 0,  A (A)  = XD(eX’l)  - L(eX*I) 


ru  in  fu 

= XH  - X eAWdy (0)  - eAOdn(0) 

- *-y*  J — r 


The  fundamental  matrix  solution  X(t)  of  (1.2)  is  defined 
to  be  the  n x n matrix  solution  of  the  equation 


(1.4) 


D (X. ) = I + f L(X  )ds,  t > 
t jn  s 


xo(0) 


Jo  , -r  < 6 < 0 


H-1, 


6 = 0 


The  following  results  of  Henry  [11]  will  be  fundamental 
to  our  investigation. 


Lemma  1.1.  If  Re  A < 6 for  all  A satisfying  (1.3) , 
then,  for  any  e > 0,  there  is  a K = K(e)  such  that 


(1.5) 


| T (t) | , |x(t)|,  |X(t) | < Ke(6+e)t,  a.e.  for  t > 0. 


Definition  1.1.  The  operator  D is  said  to  be  stable  if  there 
is  a v > 0 such  that  all  roots  of  the  equation 


(1.6) 


det  D (el)  = 0 


satisfy  Re  A < -v. 

From  the  results  of  Cruz  and  Hale  [4]  and  Henry  [11],  an 
operator  D is  stable  if  and  only  if  the  zero  solution  of  the 
functional  equation 


(1.7) 


D(y.  ) =0,  t > 0 


-IS 

zrew 
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is  uniformly  asymptotically  stable;  that  is,  there  are  constants 
K,  a > 0 such  that 

(1.8)  |yt(<j>)  | < Ke~at|<j)|,  t > 0,  ({>  e c,  D4>  = 0. 

If  D(J>  = H<J>  (0)  - J<{>  (-r) , then  D is  stable  if  the  roots 
of  the  polynomial  equation 

det  [H-pJ]  = 0 

satisfy  |pj  < 1. 

An  important  property  of  equation  (1.2)  when  D is  stable 
is  the  following  (see  [5]):  If  D is  stable,  there  is  a 
constant  aD  < 0 such  that  for  any  a > aD,  there  are  only  a 
finite  number  of  roots  X of  (1.3)  with  Re  X > a. 

If  F,G:  R -*■  Rn  are  continuous,  a nonhomogeneous  linear 
NFDE  is  defined  as 

(1.9)  — !D(xt)  - G(t)J  = L(xt)  + F (t)  . 

A solution  through  <J>  at  t = 0 is  defined  as  before 
and  is  known  to  exist  on  [a-r,°°) . 

The  variation  of  constants  formula  for  (1.9)  (see  [8]) 
states  that  the  solution  of  (1.9)  though  (a , <J> ) is  given  by 


(1.10) 


x(t) 


T(t-a)  <j>  (0)  + X(t-s)F(s)ds 
' a 

ft+ 

[dsX(t-s)][G(s>  - G(o)3 
a 


for  t > 0,  where  X is  the  fundamental  matrix  solution. 

Another  convenient  equivalent  form  for  equation  (1.10), 
is  the  following: 


(1.11) 


x(t)  - X (0) G (t)  = T(t-a)  4)  (0)  - X(t-a)G  (a) 

rt 


-r 


X(t-s)F(s)ds  - [d  X (t-s)  ] G (s)  , t > 0. 


Let  us  make  a few  other  observations  on  the  variation  of 
constants  formula  which  suggest-,  changes  of  variables  which  will 
be  useful  in  later  sections.  Let  PC  be  the  space  of  functions 
taking  [-r,0]  into  Rn  which  are  uniformly  continuous  on 
I-r,0)  and  may  be  discontinuous  at  zero.  With  the  matrix  Xq 
as  defined  before,  it  is  clear  that 


PC  = C + <XQ) 

where  <XQ)  = span  {Xq};  that  is,  any  iji  e PC  is  given  as 

¥ = * + XQb 

where  (f>eC,beRn.  We  make  PC  a normed  vector  space  by 

defining  the  norm  |i{i|  = sup  |^(8).|. 

-r<6<0 
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Let  us  define  xt(ij/)  = T(t)ij)  where  \p  e PC  and  x(iJO 
is  the  solution  of  (1.2)  through  ty.  The  operator  T(t):  PC  *> 
(functions  on  [-r,0])  is  linear,  but  T(t)  does  not  take 
PC  PC.  The  operator  T(t)  is  an  extension  of  the  original 
semigroup  T(t)  on  C.  If  we  use  this  notation,  then  the 
variation  of  constants  formulas  (1.10),  (1.11)  can  be  written 


Cl.  12)  xt  = T(t-a)  <j»  + j T (t-s)  XqF  (s) ds 


•t 

[d  T(t-s)XQ]  [G (s)  - 6(a)] 
n 59  u 


(1.13) 


xt  - XqG  (t)  = T(t-a)  [<})-X0G(a))  + T(t-s)X  F(s)ds 

■'  a 


-f 


[dgT (t— s) Xq] G (s) 


for  t > a,  $ e C.  As  usual  in  the  theory  of  functional  differ- 
ential equations,  these  integrals  are  in  Rn;  that  is,  each 
integral  is  evaluated  at  each  8 e [-r,0]  as  an  integral  in  Rn. 
Formula  (1.13)  certainaly  suggests  the  change  of  variables 


(1.14) 


xt  - XqG  (t)  = zt,  <f>  - XqG  (a)  = <|> 


from  C PC.  If  this  is  done,  equation  (1.13)  becomes 


.'.ill. ; I, I 
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(1.15)  2fc  = T(t -0)i}>  + | T(t-s)XQF(s}ds  - j [dsT(t-s)X0]G(s) 


a formula  much  simpler  than  either  (1.12)  or  (1.13).  This 
remark  will  play  an  important  role  in  the  subse.ient  discussion. 

2.  Stability  in  Nonlinear  Equations 

In  this  section,  we  give  some  elementary  results  on  the. 
stability  of  nonlinear  equations  in  order  to  show  the  previous 
transformation  from  C to  PC  can  be  of  assistance.  To  keep 
the  notation  at  a minimum,  the  most  general  results  are  not 
given. 

Suppose  F:  C-*Rn,  G:  C Rn  are  given  continuous 

functions  and  G((j>)  depends  only  upon  values  of  <j>(e)  for 

6 < 0;  that  is,  for  any  a e Rn  and  any  sequence  $n  e C, 

<{)n(0)  = a,  n = 1,2,...,  which  converges  to  <j>  uniformly  on 

compact  subsets  of  l-r,0),  the  limit  of  G($n)  exists  as 

n -»■  ® and  lim  G(<$>n)  = G ($) . The  relation 
n-H»  J 


12, 1) 


—•  !D(xt)  - G(xt)]  = L(xt)  + F(xt) 


with  D,L  as  in  Section  1 defines  a neutral  functional  differential 
equation.  Existence  of  solutions  for  initial  data  in  C follows 
from  [6] , [9] . 

The  variation  of  constants  formula  (1.13)  for  this  equation 
(2.1)  is 


wiiMiriiiiii»*riwi^M-»iniT*rii«wii>r.f'<itirii-«i)i^-iia  mwr'waar^l  C2a>,*MM-'S!,.'S 


(2.2) 


xt  - XQG(xt)  = T(t)  [<?-X0G(6)3 

+ | T(t-s)X0F(xs)ds  - |tIdsT(t-s)X0]G(xs) 


Consider  the  map 


h:  C -*•  C + <XQ> 


l 


h($)  (6)  = 


$(6) 


9 < 0 


CO)  - H ~G($),  9 = 0 


Since  G($)  does  not  depend  on  <*>  (0) , the  mapping  h has  a 
continuous  inverse,  that  is,  h is  a homeomorpnism,  (8] 

If  = h(xt),  ip  e h($) , equation  (2.2)  becomes 

(2.3)  zt  = T (t)  \p  + J T(t-s)X0F(h'1(zs))ds 


[dsT(t-s)X0]G(h  1(zs)) 


We  are  now  in  a position  to  prove  the  following  theorem. 


Theorem  2.1.  Suppose  F(0)  = 0-  G(0)  = 0 and  the  first 
derivatives  DF($),  DG($)  are  continuous  and  vanish  at  $ = 0. 


If  the  linear  equation  (1.2)  is  uniformly  asymptotically  stable, 
then  the  system  (2.1)  is  uniformly  exponentially  asymptotically 
stable. 


Proof:  From  the  hypothesis  on  G,  the  mapping  h:  C C + (XQ) 

is  a homeomorphism  in  a fixed  neighborhood  of  <J>  = 0 e C, 

= 0 e C + (Xq)  . Furthermore,  there  are  constants  k^,  k2  > 0 
such  that  in  this  neighborhood  ij>  = h(<f>)  implies  |^|  < k,  |<j>|, 
Ul  < 1^2  i ^ I * APPlY^n9  Lemma  1.1  and  the  hypotheses  on  F,G, 
we  have  there  a 6 > 0 (as  small  as  desired)  such  that  z 
in  (2.3)  satisfies 


< Ke~at  | [p  | 


+ 


tKe~a(t“s) 6k2lzs!ds 
0 


as  long  as  |z  | < e(6).  Applying  Gronwall's  inequality  to 
j zt  | ectt , we  obtain 


1 - (a-Kk9<S)  t - (a-Kk~§) t 

lxti  s lzt I ^ Ke  1*1  < Kkie  1*1 

2 

as  long  as  |zt|  < e(<$).  Since  this  clearly  can  be  assured  for 
all  t > 0 if  |<j)j  is  sufficiently  small,  we  obtain  the  result 
stated  in  the  theorem. 

As  one  sees  from  the  above  proof,  the  transformation  from 
C * C + (XQ)  reduces  the  discussion  to  an  argument  very  similar 
to  the  '-ne  for  ordinary  differential  equations.  One  could  easily 
generalize  the  above  results  to  obtain  the  more  general  stability 


0^' 
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properties  given  by  Hale  and  Ize  [10] . 

Another  interesting  remark  about  the  above  transformation 
concerns  the  manner  in  which  difference  equations  are  included 
in  neutral  equations.  Suppose  F = 0,  L = 0 in  (2.1).  Then  the 
equation  is  equivalent  to  the  functional  equation  (no  derivatives) 


(2.4) 


D(xt)  - G(x't)  = D (4»)  - G(<j>),  t > 0 


with  Xq  = (f>.  The  usual  difference  equations  are  homogeneous. 


(2.5) 


D(xt)  - G(xt)  =0,  t > 0,  xQ  = <j>. 


That  is,  the  initial  function  <j>  satisfies 


(2,6) 


D(«j>)  - G (<!>)  = 0. 


1 


■4 


It  is  a well  known  fact  that  if  the  zero  solution  of  the 
linear  homogeneous  equation 


(2,7) 


D(xt)  = 0 


is  uniformly  asymptotically  stable  (i.e.  D is  stable)  then  the 
nonhomogeneous  equation  (2.5)  satisfies  the  same  properties  if  G 
satisfies  the  conditions  of  Theorem  2.1.  On  the  other  hand,  - 
Theorem  2,1  does  not  imply  the  result  since  the  homogeneous 
linear  equation 


is  not  uniformly  asymptotically  stable  as  a consequence  of  the 
fact  D(<j>)  0 for  all  <j>  e C.  The  question  is  the  following: 

How  can  we  use  the  theory  of  neutral  equations  to  obtain  the 
above  stability  theorem  for  the  functional  equation  (2.5)? 

Let  us  now  show  how  the  transformation  h:  C -*■  C + (X^) 
solves  this  problem.  Let 

(2.8)  PCQ  = {ij>  e C + <X0)  : = 0}. 

Using  Laplace  transform  and  the  same  type  of  arguments  as  in 
Henry  [11] , one  can  prove  the  following 

Lemma  2.1.  If  D is  stable,  then  there  are  positive  constants 
K,a  such  that  if  x(ijO  is  the  solution  of 

D(xt)  = 0,  t > 0,  Xq  = if;  e C + (X0) 

then 

|x(i|))  (t)  | < Ke  at|ifi|  t > 0. 


Now  consider  the  equation 


The  variation  of  constants  formula  implies 


Xt  - XflG(xt)  = T (t)  [iJj-XqG  (<|j  ) ] - pidgTU-sjX^GUg) 


If  we  let  = h(xt),  ip  = h (<J>),  then,  a direct  evaluation  yields 


D(t p)  = D[(f>-XnG(<J>)J  = D(cp)  - G (<p)  . 


•1 


If  we  assume  <p  satisfies  (2.6),  then  D(i|>)  = 0 and 
| T (t) ip  | < K|^je  by  Lemma  2.1.  The  kernel  in  the  integral 
above  also  has  an  exponential  bound  of  the  same  type  since 
T (t) Xq  satisfies  D(T(t)XQ)  = I and  we  are  only  interested 
in  the  variation  of  T(t)Xg.  Consequently,  the  stability 
results  for  equation  (2,5)  is  easily  obtained  exactly  as  in  the 
proof  of  Theorem  2.1  after  we  have  made  the  above  elementary 
observations  about  PCQ. 

These  observations  about  functional  equations  are  the 
motivation  for  the  discussion  of  the  mixed  differential  and 
difference  equations  of  the  next  section. 


3.  A Special  Equation 


In  this  section,  we  consider  the  system 


a)  x(t)  = Ax  (t)  + By(t-r) 

b)  y (t)  - E'x(t)  - Jy(t-r)  = 
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where  x,y  are  k,m  vectors,  respectively,  all  matrices  are 
constants,  E1  is  the  transpose  of  E.  For  any  a e Rk,  <j>  e C, 
one  can  define  a solution  of  (3.1)  with  initial  value  x(0)  = a, 
yn  = <p.  If  we  define  C = C(  [-r,0]  ,Rm) , 


(3.2) 


D = 


lD2J 


: R1  x C -*■  Rk*m 


k i 


D^a,#)  = a 

D2  (a, <f>)  = <|>(0)  - E'a  - J<K-r) 


and 


(3,3) 


L = 


,0  , 


: Rk  x C -»•  RkXm 


L^(  a,<j>)  = Aa  + B<j>(-r) 


then  equation  (3.1)  is  a special  case  of  the  NFDE 


(3.4) 


D(x(t)  ,yt)  = L(x(t),yt) 


and  one  obtains  the  equation  (3.1)  by  requiring  that 


(3.5) 


^2  (**f  = 0. 


Observe  that  D(a,<}>)  = H<f>  (0)  - M<j>  (-r)  where 


14 


m 


s 

a 

H 
1 

"0  0 

II 

, M = 

<1* 

\ T J 

-E*  I_ 

0 J 

<i>  = 


Hence,  if  we  assume  the  eigenvalues  of  the  matrix  J have 
modulli  less  than  1,  then  D is  stable. 

The  characteristic  equation  of  (3.1)  is 


det  A(A)  = 0,  A(A)  = 


AI-A  -Be 


-Ar 


-E*  I-Je 


-Ar 


Equation  (3.4)  generates  a semigroup  T(t)  on  RK  x c.  If 
we  define 


(Rk  x C)Q  = { (a,<j>)  e Rk  x C:  D2(a,<j>)  = 0} 


then  (Rk  x C) q can  be  considered  as  a Banach  space.  Further- 

more,  for  any  (a,<|))  e (R  x C)Q,  the  solution  of  (3.4)  though 

k 

(a , 4> ) will  be  in  (R  x C)^  since  it  corresponds  to  the 
solution  of  (3.1)  through  (a,<j>).  Consequently, 


def 

Tq  (t)  = T(t) 


(RkxC)  Q 


(Rk  x C)Q  -»■  (Rk  x C)Q 


and  is  a strongly  continuous  semigroup.  The  infinitesimal 

Jr 

generator  of  TQ  (t)  is  j&Q  = _q^|  (R  x c)  q where 

is  the  infinitesimal  generator  of  T(t).  One  can  easily  show  that 


1 


(3.6) 


ct(_q/q)  = {X  e C:  det  A(X)  = 0} 


The  fundamental  matrix  solution  X(t)  of  (3.4)  is  defined 


I 0 
E'  I 


Xn(0)  = = H-1,  0 = 0,  Xn  (0)  =0,  -r  < 0 < 0. 

U T-l  • I*  U “ 


X11  X12 


X21  X22 


where  X^  is  a k x k matrix,  etc,  then  X must  be  a solu- 
tion of  (3.4)  with  the  initial  data  specified  above.  Therefore, 
the  matrices  X^j  must  automatically  satisfy 


(3.7) 

D2(xU(t),x: 

(3.8) 

°2  (X12(t)  'X; 

Notice  that  (3.7)  implies  xn'X2i  are  solutions  of  (3.1). 
The  functions  xi2fX22  not  sat-,'"FY  t^ie  equation  (3.1b),  but 
a nonhomogeneous  version  of  it.  However,  it  is  important  to 
notice  that  if  these  functions  were  differentiable,  the  derivatives 


would  satisfy  (3. lb).  This  is  an  important  remark  since  it 
essentially  implies  that  the  variation  of  X(t)  satisfies  (3.1). 
Now  we  consider  the  nonhomogeneous  system 

(3.9)  x (t)  = Ax(t)  + by(t-r)  + f(t) 

y (t)  - E'x(t)  - Jy(t-r)  - g(t)  = 0 

k.  in 

where  f,g  are  continuous  functions  from  10,«>)  to  R ,R  , 
respectively.  If 


’ 0 ' 

_kxm  t, 

’ f ' 

G = 

. 9 . 

e R , F = 

. o . 

then  the  equation  (3.9)  is  a special  case  of  the  NFDE 

(3,1°)  |g-  [D(x(t),yt)  - G(t)]  = L(x(t),yt)  + F(t) 

with  D,L  defined  as  in  (3.2),  (3.3)  respectively.  One  obtains 
the  equation  (3.9)  from  (3.10)  by  requring  that 

(3.11)  D2(a,«j>)  = g(0) 

If  we  let  wfc  = col (x (t)  ,yfc) , = col(a,4>),  then  the 

general  solution  of  (3.10)  is  given  by  the  variation  of  constants 
formula 


Wt  - XQG(t)  = T(t)ij>  - XtG(0) 


,F(s)ds  - [d  X(t-s)]G(s) 

' n 


As  in  Section  2,  we  extend  the  definition  of  T(t)  to 
x C + (Xq)  Y.  Then  (3.12)  can  be  written  as 

(3.13)  Wt  - XqG (t)  = T(t) [^-X0G(0)]  + | T(t-s)XQF(s)ds 

- J £dsT(t-s)XQ]G(s). 


An  element  of  the  space  Y can  be  represented  as  (a,<j>) , 

__  km 

$ = ($,b)  where,  aeR,beR,<j>eC.  Let 


Y0  = { (a,<j>)  e Y:  D2(a,4>)  = 0}, 


The  analogue  of  Lemma  2.1  for  this  situation  is  obtained  by 
using  Laplace  transform  and  applying  the  arguments  in  Henry 
[11]  and  is  stated  precisely  as 


Lemma  3.1.  If  the  eigenvalues  of  J have  modulii  less  than 
one  and  all  roots  of  (3.6)  satisfy  Re  A < -5  < 0,  then  there 
are  positive  constants  K,<*  such  that 


|x,,(t)|,  |Xp-i  (t)  | , | X. . (t)  | < Ke~at,  a.e.  t>0 


i •hi.in  iiimiilu'l  Inn iin'miiiliM'iUft 
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(3.14)  satisfies 

(3.16)  x(t)  = Ax(t)  + By(t-r)  + f(x(t)/yt) 

y(t)  = E‘x(t)  + Jy(t-r)  + f(x(t),yfc). 

Theorem  3,1.  If  D is  stable  and  there  is  a 6 > 0 such  that 
all  solutions  of  (3.6)  satisfy  Re  X < -6  then  the  zero  solution 
of  (3.16)  is  uniformly  asymptotically  stable. 


Proof : The  variation  of  constants  formula  (3.13)  is 


wfc  - XQG(wt)  = T (t)  [i|/—X0G  (i|/)  ) + 


•t 

T(t-s)XnF  (w  )ds 

J0  us 


[dsT(t-s)X03G(ws) 


where,  = col  (x  (t)  ,yfc) , ip  = col  (a ,<j>).  If  we  let 


zt  = wt  " X0G(wt}'  5 = * “ X0G(*> 

jcxm 

then  we  have  defined  the  transformation  h:R  * C -*■  Y, 
hOj<)  = 5,  which  is  a homeomorphism.  Hence,  there  are  constants 
kx,k2  > 0 such  that  |C|  < 3c1 1 ip  | , |<H  < k2  J ^ | for  |g|,|ty| 

sufficiently  small.  Also, 


DU)  = D(if>)  - G(V»)  = 0 


since  equation  (3.16)  is  satisfied.  Since  F = col(f,0),  only 


the  col(X^1/X21)  of  X is  used  to  evaluate  T(t-s)XQF(ws) . 
Therefore , we  may  apply  Lemma  3.1  and  the  same  arguments  as 
in  the  proof  of  Theorem  2.1  to  complete  the  proof  of  the  theorem. 

It  is  obvious  that  one  could  generalize  the  results  in 
this  section  in  many  ways.  The  perturbations  g,f  could  depend 
upon  t as  long  as  all  estimates  are  uniform  in  t.  Also, 
and  more  importantly,  the  linear  equation  can  be  much  more 
general.  In  fact,  we  could  have  considered  an  equation  of  the 
form 


where  L = 
functional  and 


R 


~ D(x(t),yt)  = L(x(t),yt) 

x c R^xm  is  an  arbitrary  continuous  linear 


D(a,<}>)  = H 


a 

4>(0), 


- dx($) 


where  H is  a (kxm)  x (kxm)  nonsingular  matrix  and 


r° 

dx(<j>)  = J [dy(O)]<J>(0) 


with  y an  m x m matrix  of  bounded  variation  with  y(0)  = y(0 
This  means  the  results  apply  to  the  more  general  equation 
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x(t)  = L1(x(t)/yt)  + f(x(t),yt) 

D2(x(t)/yt)  = g(x(t),yt) 

where 

D(a,<J>)  = col(a,D2  (a, 4>) ) . 

These  same  remarks  apply  equally  as  well  to  the  results 
in  latter  sections. 

4 . Absolute  Stability 

Let  D,L  be  as  in  section  3;  let  h > 0 be  given  and 
let  £:  R -*■  R be  a given  continuous  function  satisfying 

(4.1)  h^cr2  < af(a)  < h2a2,  0 < h^  < h2  < h. 

Jc  m 

Let  c be  a k-dimensional  row  vector,  b^  e R , b2  e R 
be  constant  vectors  and  consider  the  system 

(4.2)  x(t)  = LJL(x(t),yt)  + b1f(c(t)) 

D2(x(t),yt)  = b2f(a(t)) 

a = cx. 

Our  objective  is  to  apply  the  method  of  Popov  to  determine 
sufficient  conditions  for  the  absolute  stability  of  system  (4.2). 
The  variation  of  constants  formula  implies  the  solution  of  (4.2) 


through  ip  = (a,<p)  e R * C satisfies  the  equation 


(4.3)  wt  - XqG  (a(t) ) = T(t)  [^-XqG(ct(0))]  + T(t-s)XQF{a(s))ds 


ft 

))  + M 


-r 


. (dT(t-s)XJG(o(s)) 
J n “ u 


where  = col  (x(t)  ,yt) , G = col(0,b2f),  F = colfb^f^O)  and 
a = cx. 

As  in  the  previous  sections,  let  us  simplify  the  formula 

(4.3)  by  putting 


(4.4)  zfc  = w - X0G(cx(t}). 

If  we  let  = (u(t),vt)  then  (4.4)  is  equivalent  to 

(4.5)  u(t)  = x(t) 

vt  = yt  " X22,0  b2f(cx{t)) 

or 

(4.6)  u(t)  = x(t) 

vt(0)  = yt(0),  0 < C 

vt(0)  = yc(0)  - b2f(cx(t)) 


Equation  (4.3)  for  z.  becomes 


(4.7) 


zt  = T(t)z( 


■s)XQF(a(s)  )ds 


■s) XQ]  G (a  (s) ) 


since  <j(t)  = cx(t)  = cu(t). 

For  the  following  discussion,  it  is  convenient  to  have  this 
equation  (4.7)  explicitly  in  RkXm.  It  is  easily  seen  that  (4.7) 
is  equivalent  to  the  equations 


(4.8) 


a)  u(t)  = u (t)  + j (t-sjb^f  (a (s)  )ds  - 

“ j^WsX12(t~s)3b2f(cr(s)) 

b)  v(t)  = v°(t)  + | X21  (t-sjb^  (u  (s)  )ds 

- j [dsX22 (t-s) ]b2f (c(s) ) 


By  Lemma  3.1,  we  know  that  if  D is  stable  and  there  is 
6 > 0 such  that  Re  X < -<$  for  all  X e <?(J3^)  then  T(t)Zg/ 
that  is,  u°(t),v°(t),  and  xi;L  (t)  ,X21  (t)  ,3^  (t) , i,j  = 1,2 
approach  zero  exponentially  as  t -*■  «■ . 

Consider  (4.8a).  As  a(t)  = cx(t)  = cu(t),  we  have 


cr(t)  = cu°(t)  + | cX^  (t-s)b1f  (ct(s)  )ds 


-t  . 

O',,  (t-s)b,f (c(s) )ds. 
n 1 


Let  cu  (t)  = y (t) , cX11(t)b1  = k.^t),  cX19^t)b9  = k,  9 (t) 


12  v '2  12 


H jt*  *iw«  -v  mm  , 


"’9>  alt)  = Y(t)  + J0"kll(t-=)f  («(•))  as  + }\2(t-s)f(a(s))ds 


Notice  that  v(t),  lliil  v /+-\  t ^ 

r dt  ' kn  / ^22  ^ tend  to  zero 

exponentially  as  t oo. 

Now,  define 


fT(t)  = 


f (cr  (t) ) / 0 < t < T 


t > T 


cfT(t)  = ^ 


a (t)  / J < t < T 

o / t ;>  t 


YTCt)  - 


Y (t)  / 0 < t < T 


i ; 


fT  rT 

J0kll(t"S)f (0(8)) as  + j k12(t-S)f (0(s))dS/  t > t 


Thus 


4.10)  <;T(t)  - YT(t)  - |0k1i(t-s)fT(a)ds  + f k12(t-s)f  (g)ds. 


T f 

. (4.11)  X (T)  = Jcr(t)  - y(t)  -If(a(t)) 

JOt 

= . „■  °t ' ‘ h fT  + q(d?  - a?]  > fTdt 
« ^ 

.00 

3.1  N fw  1 ^ 

“ 2 7f  I E c^qi  T h ij>  ~ Y«p)  } fipdto 

* *-*00  ~ 

where  N denotes  the  Fourier  transform  and  f*  is  the  transpose 
of  f.  From  (4.10),  by  the  convolution  theorem,  we  have 


^ll’^T  + *12*T  *lA  + iw*12*T  * 


Hence,  (4.11)  becomes 


if”  ~ » 

X(T)  = 7*  {Re[(l+iwq)(ki;L+ittk12)]  - £}|f  pdw. 
— 00 


If  the  condition  Re  L (l+itoq)  (k1i+^a)kx2^  ~ £ < 0 is  fulfilled, 
then  x(T)  < 0*  This  means  that 


£ f (a  (t) ) + q SgMl  f(c(t))dt 


< |^|y(t)  + g f(^(t))<3t- 


lim  f (tf (t)  = 0.  This  implies,  by  (4.8a),  that 

t->-<n 

lim  u(t)  = lim  x(t)  = 0.  From  (4.8b)  it  follows  that 

t-H»  t^00 

lim  v(t)  = 0,  since  lim  f (a (t) ) = 0.  Hence  lim  y(t)  = 0. 
t~>-<»  t"^00  t^00 

Thus,  the  following  theorem  is  proved; 

Theorem  4.1.  If  D is  stable,  all  roots  of  (3.6)  satisfy 
Re  \ < -<S  < 0 and  there  is  a q > 0 such  that 

Re[(l+iwq)  (k^+iak.^)  ] ~ E - 0 

where  k.^  = cX^b^,  k12  = cXi2b2'  then  the  system  (4.2)  is 
absolutely  stable  for  every  function  f satisfying  (4.1). 

Rasvan  [12]  has  studied  the  same  problem.  He  gives  the 
condition  on  q in  terms  of  the  transfer  function  n (iw) . 

If  we  compute  + iwR.^  we  see  t^iat 

n(i«)  - k^^  + iwk12  ' 

One  could  easily  study  the  first  critical  case  using  the 
same  arguments  as  in  Halanay  [7]  and  Datko  [3] . 
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